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^■-^ \ Abstract 

We show that (1 + 2) nonlinear Klein-Gordon equation with negative coupling 
admits an exact solution which appears to be the linear superposition of the plane 
wave and the nonsingular rational soliton. We show that the same approach allows 
pi ^ to construct the solution of similar properties for the Euclidean (\)^ model with 

D \ broken symmetry. Interestingly, this regular solution will be of instanton type only 

in the D < 5 Euclidean space. 



> 

^ . 

^ '• 1 Introduction 

in 

^ ' The nonlinear Klein-Gordon equation (0 - or |0| -models) belongs to the class of those 

Q ■ equations that frequently arises in seemingly different physical applications: see, for ex- 

C^T ■ ample, [1] (cosmic strings), [2] (field theory), [3J (condensed matter physics, e.g. liquid 

crystals). 

Unfortunately, this important equation is not an integrable one. In the particular case 
^ ' of (1 + 1) one can still obtain stationary solutions (for example, kinks for the 0^-model), 

■ but in the case of general position (1 + D) there are no general methods allowing one to 

construct the exact solutions. There are some exceptions, however: most notably in the 
cases when the required solutions have sufficient number of symmetries. Such symmetries 
often allow to reduce the problem to the one- dimensional one and to successfully integrate 
the equation. For example, in the massless theory with the potential 

VW = -'4. (1) 



and upon introduction of the special boundary conditions the Euclidean solution with 
0(4) symmetry can be obtained. This solution, the so-called Fubini instanton is one- 
parametric, nonsingular, and has the form 0(r) (r = -y/aJf, the Xj are the Euclidean 
coordinates) and finite Euclidean action Jl]. A somewhat surprising fact is that for any 
non-vanishing m and the potential 

= (2) 



there are no instanton solutions in the theory. This is the good example of the nontrivial 
properties of solutions for nonlinear multidimensional equations like the mentioned 0^- 
model. 

In this article we are going to construct the generalization of the Fubini instanton 
(FI), i.e. we are going to find out the solution of the Euclidian equation of arbitrary 
dimensions and with the non-vanishing m. As we shall see, this task is accomplishable 
for the generalized potentials of the kind: 

V{^) = Vo + a<P + + - 

To be instanton this solution should be generated by the very special boundary conditions 
and must result in finite Euclidean action. As we shall see this will be the case only for 
the Euclidean spaces with D < 5. 

In order to avoid any possible confusion, let us make one little remark: everywhere 
throughout this article the A (coupling constant) is assumed to be a positive quantity 
(A > 0). Since the models we are going to discuss will (with just a couple of exceptions) 
all have the negative coupling constants, our restriction will just mean the presence of 
the negative sign in front of the corresponding terms (see, for example, equations ([T]) and 

(ED)- 

The plan of the article is as follows: first of all, we are going to consider the {1 + D) 
1 01^- model in the Minkowski space (summation is implied over repeating contravariant 
and covariant indices): 

a^a'^0 + m20- A|0p0 = 0, (3) 

with the metric 

g^^ = diag(+l, -1, -1, -1). 

As we shall see, in the case D = 2 the equation ([2D admits the nonsingular solutions 
4>{x^), /X = 0, 1, 2 such that 

|0(x'')| B = const at + oo. 

We assume that similar solutions can be obtained for the case D > 2 also, but this is still 
a hypothesis. 

Of course, at first sight the models with negative coupling seems to bear little to 
no connection with the reality. In fact, since the potential is negative, the action is 
not bounded from below and therefore, it might not be stable quantum mechanically. 
However, it appears that models with the negative potentials can be very important in 
such areas as cosmology and strings theory [5]. Thus, the problem to find interesting and 
nontrivial exact solutions of such models should be regarded as an actual one. 

We are discussing such solutions because they somehow hints the correct anzats for 
the model we are interested in with the real scalar filed in the euclidean space. We will 
show (in Sec. 4) that by the usage of the very simple ideas it is possible to find out the 
generalization of the FI - which is specifically the main aim of this article. This solution 
is valid for the Euclidian space of arbitrary number of dimensions and arbitrary sign in 
front of A. However, to obtain a regular (i.e. without any peculiarities) solution one need 
to choose the negative sign. 

The most amazing fact here is that of all the possible dimensions of the model, the 
dimension D = 4 stands out. In particular, case D = 4 is the only one that requires a 
special treatment, which treatment we are going to discuss in section 5. 



2 Equations 



In [6] the novel exact solution of the Davey-Stewartson II (DS-II) equations describing 
the soliton on the plane wave background is obtained. This solution has the form 



(j){x,y 



t) = Se*'(^'^'*) 1^-1 



(4) 



where s{x, y, t) and Pi (a;, y, t) are linear functions whereas P2{x, y, t) is polynomial of the 
second order such that P2{x,y,t) > for any values of x, y, t. 

The aim of this work is to show that the nonintegrable (1+2) |0|^ with negative 
coupling admits the similar solution (see (fT^ ) with 



Pi = a/^x^ + a, P2 = ij^^x^'x" + + A^, 



SnX^ 



where 



{A^y = A^, 



B* = P, 



The solution fll4p will be nonsingular if 



V12 V22 



> 0, 



^700 ^01 ^702 
Vol Vn V12 
V02 r]i2 V22 



> 0, 



bo hi 62 

&o ^700 Vm V02 

bi r]oi Vu V12 

h V02 rji2 V22 



> 0. 



Substituting ([5]) into ([3]) one gets: 

^P2 [P2 (Pi - P2) d^d^^s + 2J^^d^s] + P^d^J^^ - 2J'^d^P2+ 
+ (Pi - P2) [-AP2(Pi - P2)(P2 - Pr) + (m2 - d^sd>^s)Pi] 

where 

J^' = P^d'^Pi - Pid^Pi. 



0, 



(5) 



(6) 



(7) 



Using ([7]) for the particular case 6/^ = one will end up with the following system 
Sf.s^' + - m2 = 0, 

(2isX - ^^^(a + «*))^^ - a(-A52(a + 2a*) + 2'qt'^)A^ - \B'^\a\'^a = 0, 
«M + = 0) 

{2is^a^ — XB^{a + a*)) rjaf^ — Ais^^ri^adis — XB'^aaap = 0, 



+ etc 



0, 



- (a + 2a* - 2A^{a + a*)) XB'^ + SiaA^'q^^s^a^ Vai3 + AS^ (A^ + a* - 2a) a^ap- 
-Sarj^aVis + 4:iA'^s^'r]f,c,ap = 0, 



-XB^ {2a* - 2|a|2 + a^) + 2ril^A^] a„ + AA^ {a^' + ias>') r/^„ = 0. 

Thus for the case {1 + D) one should solve the system of 

^^^^ ^ iD + 2){2D' + 5D + 7) 

algebraic equations with two additional conditions: 

(i) the inequalities ([6]) must hold and 

(ii) level lines of the function P2 must be the closed curves. 



3 D = 2 Solutions 

In the case D = 2: N{D) = 50. Assuming tjq, a, (3, 7, p and b to be the arbitrary real 
parameters, define three Lorentzian vectors 



such that 



p2 sin[2(^-a)] 

r]i = r]o cos a, r]2 = r]o sm a, a = — -, 

4?7o sm(a; — 7) 

^i = pcos/3, ^2 = psin/?, ^0 = pcos(/9 - a), 

a cos{P — 7) 



(9) 



6*1 = (jcos 7, 02 = asm 7, 60 



cos{P — a) 



and besides 



2-2/ m 2 cr2 sin(2/5 - 7 - a) sin(7 - a) 
A = 2p sm (a — p), m = 



cos^(/5 — a) 



(10) 



Then the solution of the system ([8]) such the (EI) and two conditions (i) and (ii) holds has 
the form: 



Vt.u = - "^b^^Vu + 4 (62 + B') 7]^r],, = {2B^ - b^) r]^ + b^, + 6^ 



(11) 



'2B 



This would be the case if 6^ = 0. If 6^ 7^ then one have a shghtly more comphcated 
system them ([8]). The particular solution can be obtained, though: 

a = 4 (1 + ^(c^ct - c^co) , A' = AB\x' + ^1 ^^^^ 

i{cicl - C*C2) 
2|ci|2 



K Ci 




C1C2 + C*C2) 


25 


Cl 


2 



where Ci^2 1^ = =tl are the arbitrary complex constants and is similar to the previous 
one (see ( ITT!) ). 

In order to show that (EI) and two conditions (i) and (ii) hold it is sufficient to show 
the following: 

r/^.x'^x" + + A^ = 1 + AB^ [{-^^.x^ - xf + iK^" + 
where 

Remark 1. In articles [Z], [H] we has presented the method of construction of the 
exact solutions of equation ([3]), based on the usage of the self-similar change of variables, 
effectively reducing ^ to the nonlinear Schrodinger equation. It is possible to show that 
each solution of NLS corresponds to the solution of ([3]) (the contrary statement is certainly 
false). In particular, using the multi-soliton solutions of the NLS one can find out the 
corresponding "multi-soliton" solutions of ([3]). Geometrically, such solitons are similar to 
those of the Kadomtsev-Petviashvili equation with the streamlines being just a familiy of 
a parallel straight lines. Unlike to the solutions, constructed in the works 18] , the (JT4l) 
appears to be the linear superposition of the plane wave and the rational soliton with 
level lines having the form of ellipses on the plane xy. 



4 Generalized Fubini instantons in the 5D model with 
broken symmetry 

We have already noted that in the 0^ model with potential ([2]) for any nonvanishing m 
there exist no instanton solutions. On the other hand, we have shown that the complex 
model (l3l) does have the nonsingular solutions with the following behaviour: 

|0(a;^)| const at + y"^ ^ 00. 

Taking in regard the similarity between ([2]) and ([3]) models (both of them has the nonlin- 
earity of the fourth order) it is natural to search for the regular solutions of the real- valued 
model 0^ with the similar behaviour. In this section we will show that such solutions in- 
deed exist and that their Euclidean action is bound to be finite, i.e. such solutions are in 
fact instantons. However, this would be true for the 5D space-time only. 



The D=4 Fabini instantons has the form 



X 2-^' (13) 
A + 



where p is arbitrary. The corresponding Euchdean action is 

^2 



s 



(4) 



3A 

Such solutions are interesting in studies of the the tunnehng trajectories. All we need 
now are two boundary conditions: 

(i) 0(r = oo) = 00 and 

(ii) d(t)/dr\r=Q = 0. 

Instantons (I13p are the solutions whose minimal action posses the 0(4) symmetry of 
Euclidean space and which are created by the boundary conditions (i) and (ii). 

In the previous sections we have shown that ([3]) admits a certain type of solitons 
with the boundary condition akin to (i). Although the calculations so far has been done 
only for the (1+2) case, it seems natural to expect the similar solutions in the higher 
dimensions. Let's look for the generalization of f[T^ having the following form: 

^i^') = h + . (14) 

where the Euclidean variables has been renumbered with the aid of the Latin indexes 
= 1...D). Unfortunately, the direct check here would be extremely cumbersome, so 
we'll restrict ourselves to the case of the 0{D) spherical symmetry. In this particular case 
the Euclidean equations will have the form: 

^ D-l dcj) _ dV{(j)) ^^^^ 
dr"^ r dr dip 

Since the number of the free parameters will get too restricted, we'll add such parameters 
by introduction of the broken symmetry into the model: 

V{4>) = + a0 + ^ + ^ - (16) 

Using the ansatz f|T4|) by complete analogy with the ([3]) we'll end up with the required 
solution. 0: 

40o(^-4) 



p \ + p-^) 



where 



2_2A0|CD-4)2 _ 0o(A0g + m2) 

P 2 ,0 



p' ' 2 ' 

/3=— , p = 'i\<g^m\ 

200 

where 0o is arbitrary. One can see that both (i) and (ii) are hold. The calculation results 
in 

02(3A02 + 2m2) 




12 



^We omit most of the calculations made in this sections due to the fact that they are essentially the 
same as the previously performed ones 



so in order to obtain the finite Euclidean action one should take 



Vo — — 



12 

Then the Lagrangian will have the following asymptotic: 

2 




-6 



when r ^ oo. In order to calculate the action one has to multiply the Lagrangian by 
r^^^. It is easy to see that the Euclidean action would indeed be finite only in the case 
D < 5. It seems that D < 4 cases are physically meaningless and provide little interest 
for the quantum cosmology (see, however, the note 2). Case D = 4 requires the specific 
treatment and will be considered in the next section. Hence, the only case left for us to 
check out is the one with D = 5, whose Euclidean action has the form: 



S 




(5) _ 87r2 /•~^4 / 1 fd(f)Y , ^^^^^ \ 871^00 72 



and the probability of the bubble formation can be estimated as 



P ~ exp 



Remark 2. Multiplying the ( fT7|) by the oscillating exponents and making such choice of 
the parameters as to set the cubic term {f3 = 0) to zero, we'll end up with the solution of 
the non-homogeneous multidimensional NLS: 

i^t + A(^)^ + (k- SAc^) ^ + A|^|2^ = -2AcV''*, (18) 



where A*^^) is a D-dimensional Laplacian and k, c are arbitrary constants. An exact 
solution of such equation will be of a form: 




By a particular choice k, = 3Ac^ it is possible to reduce the left hand side of (|T8l) to the 
traditional NLS form. This solution works for both focusing and defocusing NLS, but for 
the latter the equation will have a singularity. It is well known that in case of D > 1 
and under the assumption of the non- vanishing right hand side the equation f|T8|) becomes 
an non-integrable one, posing quite a challenge for those who dares to find out its exact 
solutions. On the other hand, NLS is the one of the most commonly used equations in 
physics. So we hope that the existence of the exact solution (lT9l) for equation (fT8l) would 
be useful in different physical applications. 



5 Case D = A 

The Bogolyubov transformation — + c is the one commonly applied for the models 
with the spontaneously broken symmetry, i.e. for the case when either Lagrangian or 



Hamiltonian allow for some kind of symmetry, while the ground state is not invariant 
with respect to the corresponding symmetric transformation. However, the Bogolyubov 
transformation can also be applied in a slightly more general context. Lets take for 
example the solution (1171) . Since as r — oo, — (pQ, it is quite natural to try a new field 
variable $ = — 0O; which goes to zero as r ^ cxd. 
Lets consider the following simple potential 

= ^ + K,', (20) 

where lo'=const (we remind our readers, that A > 0, i.e. the model under description is 
the one with the positive coupling constant). Introducing the new variable 0, we end up 
with 

V{<P) = ^ - c<Po<P' + ^<P' - A0^$ + V^. (21) 

Comparing (12T|) and (fT6l) (with regard to the A —A substitution), we deduce the 
expression for the mass 

m = V3A0O- (22) 

The direct substitution of fl22l) into f|T7|) will obviously result in singularity. In order to 
avoid this, let us point out that for D = 4 the expression contains the indeterminacy, 
whose elimination might in principle give us the reasonable result. Hence, the "symmetry 
breaking" procedure □ is possible for the D = A only. 
Now introduce the following handy substitution: 

- 3Ac0^ = e{D - 4) (23) 

and take the limit e — * 0. Substituting (123!) into (fT7|) will result in: 

^ = <Po = -rPhr.- (24) 

^2^2 — 2A0O 

Keeping e as a finite quantity, it is possible to substitute fl23l) directly into f|T5l) with the 
potential fl?Il) and to convince oneself that for D = 4 the fl25]) is the solution indeed. Let 
us emphasize that this solution corresponds to the case of nonzero mass, but is singular 
and, hence, leads to certain problems with calculations on the Euclidean action and the 
probability of the bubble formation. However, there are still some ways to bypass those 
difficulties. First, if we'll assume in ( 121]) Vq = A0q/4 it will give us 



0. 




Second, calculating the lagrangian 



320ge^ {rh^ + 2Ac^) 
(r2e2 - 2Ac2) 



^This can be understood literally: while initial model ([20]) has had the discrete reflection symmetry 
<I> — <I>, the pT|) does not. However, it is more convenient to refer to this as to the "hidden", not 
"broken" symmetry: after all, the later term traditionally has a precise and strict definition - see the 
beginning of this chapter. 



and substituting it into the expression for the Euchdean action 



= 2t:^ r r^Ldr. (25) 
Jo 

will end up in expression, having the inherent singularity (pole of the fourth order). How- 
ever, the troublesome term can be effectively separated from other ones by representing 
(ESD as: 

Sf = + Sg), (26) 

where 

Note that the probability of the bubble formation is defined by 



^vr2 /87r2 



exp 



^(^) UiV^exp pi- . (27) 



3A °^ / " \ 3A 



Note also, that the emergence of the "strange" multiplicative constants, such as A^oo and 
alike during the calculations of amplitudes by the functional integral method is quite com- 
mon and should not be highly surprising [9]. Following the generally accepted (although 
not strictly proven) methods, one can simply cancel the N^o quantity in the P. As a 
result, we'll get the expression 



exp 



which is similar to the expression, obtained for the normal Fubini instanton [TO] up to the 
sign . 

The problems, described above are the inherent ones for the models with the positive 
coupling constant. In case of the negative coupling constant A ^ — A the four-dimensional 
action has the normal form 

Sf ^ (28) 

and provides exactly the same probabilities as the massless model with the Fubini instan- 
ton. 



6 Conclusion 

There are three open questions which are interesting problems for the future investigation. 

The first question is the generalization of the method to the case D > 2. At a 
moment there are no obvious reasons to prevent us from using our approach in the multi- 
dimensional case. The only price we would have to pay working there would be the rapidly 
growing amount of calculations. In fact, the number of algebraic equations N{D) which 
should be solved grows as D^: N{3) = 100, A^(4) = 177 and A^(IO) = 1542. Our point 
about this is following: even if it can't be done by hand, the problem can still be resolved 
via the computer programs (such as MAPLE or Mathematica) , powerful enough to handle 
such type of calculations. 

The second question is the necessity of the negative sign in front of the corresponding 
terms with A for the case of |0|^ model (Sec. 2 and 3). Is it possible to construct exact 
solution with positive sign (i.e. with positive coupling)? The answer is yes, but in this 



case we'll have to deal with the singular solution. This singularity is located along the 
hyperbola and the physical meaning of such solution is unclear for us. 

The last question is connected with the amazing fact that the dimension D — A stands 
out. Whether there is a simple explanation of this fact? 

We conclude that one can use the generalized Fubini instantons in quantum cosmology 
only for the case of the single infinite extra dimension. 
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